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Abstract: 

A probabilistic  analysis  of  the  fatigue  crack  growth  for  reliability  growth  calculation  on 
the  mechanical  component  is  presented  on  the  basis  of  fracture  mechanics  and  theory  of 
random  process.  The  loading  is  postulated  to  be  stationary,  narrow-band  random 
Gaussian  process  and  consequently,  randomized  Paris-Erdogan  law  is  applicable.  As  a 
specific  problem,  a thin  plate  having  a central  crack  is  analyzed  by  two  analytical 
methods  “stochastic  averaging”  and  “damaged  linear  accumulation”.  Simultaneously, 
the  aforementioned  plate  is  being  analyzed  by  the  utilization  of  “NISAII”,  Software. 
Analysis  made  by  “NISAII”  is  on  the  basis  of  Monte-Carlo  simulation  random  analysis. 
In  the  end,  all  the  results  are  compared  with  each  other  and  conclusion  is  drawn. 
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Introduction: 

One  of  the  main  contributors  to  the  structure  failure  of  an  aircraft  is  a phenomenon 
known  as  fatigue.  Most  of  the  mechanical  components,  especially  aircraft  structure,  have 
been  designed  on  the  basis  of  “Fail  safe”  and  “damage  tolerance”  concepts.  Providing 
assurance  of  no  failure  operation  is  highly  important  in  these  components.  Reliability 
evaluation  can  provide  useful  information  in  the  fatigue  control. 

Evaluation  of  fatigue  effect  on  mechanical  components  containing  initial  crack  is  rather 
complex.  The  complicity  arising  from  crack  growth  in  component  that  depends  on 
essential  properties  of  material  under  fatigue  load  also  depends  on  geometry  of 
component  as  well  as  environmental  conditions.  In  the  laboratory  controlled  conditions. 
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if  an  experiment  is  repeated  many  times  over  and  test  conditions  remain  the  same,  the 
results  are  not  deterministic,  and  in  most  cases  are  scattered.  Critical  elements  designed 
on  the  basis  of  the  above-mentioned  points  will  be  very  reliable  One  of  the  important 
concepts  in  reliability  analysis  is  the  definition  of  failure  criteria,  which  can  estimate 
probability  of  useful  operation.  In  mechanical  components,  existence  crack  is  one  of 
major  failure  causes,  hence  reaching  crack  length  to  a critical  length  can  be  defined  as  the 
failure  criteria.  Equations  of  fracture  analysis  are  in  the  form  of  probability,  therefore  we 
can’t  apply  deterministic  fracture  mechanics.  Stochastic  behavior  should  be  evaluated  by 
statistical  and  probability  theories  for  reliability  analysis  which  introduce  complicity  in 
analysis  results.  Reliability  can  be  analyzed  with  obtained  results  and  proper  modeling 
and  solution  by  stochastic  method  Thus  becomes  an  important  point. 

In  this  regard,  there  are  a host  of  stochastic  methods.  Moreover,  recently  presented 
methods  reduce  essential  complicity  of  probability  calculations.  Stochastic  averaging 
method  for  fatigue  crack  measurement  in  a component  under  random  loading  presented 
by  Zhu  and  Lin  (1992)[1]  is  to  name  one  B.Kececioglu  in  1998  analyzed  reliability  of 
mechanical  component  under  wide-band  and  narrow-band  loading  by  definition  of 
damage  function  [2],  Results  obtained  from  the  aforementioned  researches  are  used  as 
comparative  criteria  in  the  presented  study  here 

In  this  paper,  reliability  analysis  of  a mechanical  component  is  evaluated  under  critical 
conditions.  In  addition,  failure  cause  is  assumed  to  be  the  fatigue  fracture  that  arises  from 
crack  growth.  Component  behavior  is  evaluated  on  the  basis  of  probability  fracture 
mechanics  and  random  vibration  theories.  Critical  conditions  are  presumed  to  be  simpler 
assumptions  for  reducing  governing  probability  complicity.  On  the  basis  of  these 
assumptions,  the  problem  is  designed  and  solved  by  governing  theories.  Property  of  a 
dynamic  system  is  light  damping  and  stationary  narrow-band  Gaussian  random  loading 
process.  Consequently,  the  system  response  would  be  come  close  to  the  system's 
resonance  response,  (i.e  high  stress  fatigue  is  taken  into  account ) “Stochastic  averaging 
method”  and  “damaged  linear  accumulation  method”  is  used  as  analytical  solutions. 
Software  analysis  (by  NISAII)  is  presented  on  the  basis  of  Monte-Carlo  random 
simulation  method. 

Methodology: 

In  this  section,  analytical  and  software  methods  are  used  for  the  reliability  growth 
determination  in  a cracked  mechanical  component  under  fatigue  failure.  Considering  the 
complicity  of  methods  for  reliability  growth  in  this  field,  it  becomes  necessary  to  design  a 
problem  with  simpler  conditions  in  order  to  make  it  possible  to  expand  into  general 
conditions. 

Problem  Assumption: 

Dynamic  Analysis 

• one-degree  of  freedom  and  light  damping  system 

• stationary  wide  band  Gaussian  random  loading  process(white  noise) 

Fracture  Analysis 

• only  one  crack  in  component 

• length  and  location  of  crack  known 
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• material  behavior,  brittle,  homogeneous,  isotropic 

• high  stress  fatigue 

Reliability  Analysis 

• failure  case,  fatigue  fracture 

• failure  criteria,  a=ac 

Model  Selection: 

Model  assumptions  should  be  one-degree  of  freedom,  light  damping  and  wide  band 
loading.  Test  specimen  concerned  is  thin  plate  with  a central  crack,  however,  it  can  be 
expanded  to  other  specimens. 

Problem: 

Consider  a thin  square  plate  1 x 1,  shown  in  Fig.  1,  with  an  initial  central  crack  of  length 
2ao  and  supporting  an  infinitely  rigid  heavy  mass  M at  its  end.  The  plate  is  idealized  to  be 
massless,  homogeneous,  isotropic,  and  with  light  damping,  and  the  mass  Mis  subjected 
to  a Gaussian  load  process  perpendicular  to  crack  (Y)  with  a wide-band  one-sided 
spectral  density  G. 


5(t) 


Fig.  1 Thin  plate  containing  central  crack  and  mass  in  its  end 


Analytical  Solution  Methods: 

For  calculation  of  plate  reliability  growth  presented  “stochastic  averaging”  and 
“accumulative  linear  damage”  as  analytical  methods 

(l)Stochastic  Averaging  Method[3] 

Model  of  System  Dynamic 

For  one  degree  of  freedom,  the  system  can  be  considered  as  dynamic  model  of 
mass,  spring  and  damping  system,  then  the  dynamic  equation  will  be: 

MX(t)  + CX(t)  + g[A(t)]X(t)=4(t)  (1) 

The  stiffness  function  g[A(t)]  can  be  approximated  by  following  polynomial 
expression[4] 

g(a)  = g(0)(l  -1.708«2  +3.081u4  -7.036«6  +8.928«s  -4.266a10)  (2) 

2a 

Which  g(o)  is  without  crack  plate  stiffness  and  g(ag  = 0)  = E.  A / L and  u = — . 

We  shall  assume  that  the  stress  is  equal  to  the  displacement  multiplied  by  a constant. 
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Fracture  Analysis 

Based  on  Paris-Erdogan  model,  the  equation  for  crack  growth  rate  is: 

dA  = W{A)  «(A*y 

dt  2k 

(3) 

Where: 

tsK  = 2]j;!Gg(ah(a)R 

(4) 

R , stress  range,  an  approximate  expression  for  h(a)  has  been  provided  in 

[4]as  following: 

h(a)  = V«(0.467-0. 514«  + 0.960m2 -1.421m3  +0.782i/4) 

(5) 

now  eq.(3)  can  be  rewritten  as  follows: 

~ = yQ(a)s^ 

dt 

(6) 

Y~{2j(*Glc))f 

(7) 

Q(A)  = a>Wh(A)Jte(A))]> 

(8) 

In  which  S is  the  stress  envelope  process  and  R (t)=2  S(t).  Then  by  solving  the  above 
equations,  the  following  transition  probability  density  of  crack  size  A (t)  will  be  yield 


q(a,t\a0, 0): 


yJ(27a)Q(a)cr<t>(myfl\^T) 


exp 


f du 

mt 

io0(") 


a >an 


(9) 


Where: 


m = y T(1  + p\2) 

(10) 

<r2=Y2±clzx 

n~\ 

(11) 

(~1)tw!  r(\+k  + p\2) 

(n-k)\(k\)2  1 

(12) 

=2$”  p2"(r)dr 

(13) 
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p=exp[-C  t /(2M)]  correlation  coefficient  of  loading  process. 
The  reliability  function  follows: 


<X> 


mt  - 


R(t,acr\a0,0)  = £ ' q( a, t\a0 , 0)da  = 1 — 


a -Jt 

'4t\a 


(14) 


Where: 


_ r»t,  da 
Zc,~k  Q(a) 

(2)  The  Linear  Accumulation  Hypothesis  in  random  fatigue  crack  growth 
Paris-Erdogan  model  is  given  in  the  form  [5] 

= C(Ak)m  (15) 

dn 

Where: 

Ak  = ASVrar|(a) 

When  the  crack  size  is  very  small  in  comparison  with  the  component  dimensions,  the 
above  equation  can  be  rewritten  as 

d°-  = CG(a)(tsSy  (16) 

dn 


Where: 


G(a)=[ri(a)V^]m 


and  Eq.(16)  can  be  written 


'#'(«)=/■ 

O0 


dz 

G(z) 


(17) 


Where 


dY(a) 

dn 


= C(AS)m 


the  damage  indicator  is  defined  as: 


D = 


'i’M 


(18) 
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It  can  be  shown  that  when  a=aO,  D=l;  and  when  a=ac,  D=1  The  increase  rate  of  the 
damage  indicator  therefore  is 


dD  2 mfC(ASy 
dt  ~ 

As 

Where  S is  the  stress  amplitude  process,  and  S = — . 


(19) 


Definite  the  coefficient  A as 


a l 

A = <o=;w 

2 m 

C 


(20) 


It  can  be  shown  in  [5]  that  the  damage  accumulation  under  stationary  stress  loading 
process  is  a normally  distributed  random  variable;  i.e., 


£>(f)oc  N[D(t),a Dm] 


For  narrow  band  stress  process: 


D(0  = -i-(V2OrJftr(-  + l) 

2 7t  A 2 


{4l<jx)b  6 | la>„  (p]  (b) 
A '2  y 2k  t; 


(Pl  (b)  which  is  given  in[2]  and  a and  P definite: 

<*  = -!=  ,P=~ 

Vw  u 


(21) 

(22) 

(23) 

(24) 

(25) 


(V2crJ6r(^  + l) 


V 


a>n 

2 k 4" 


The  reliability  is  then  given  by 


/f(0  = <D{-[(^),,2-(^)''2]) 
a p p 


(26) 


(27) 
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Analysis: 

Reliability-Time  curve  shown  in  Fig.2  for  analytical  method. 


Fig.2  curve  of  Reliability  -Time  with  two  methods 


Physical  Data  of  Problem: 

2ao=0.00254  [m] 

2ac=0. 01016  [m] 
C=3502.536  [Kg/S] 
E=68.95  [GPa] 


G=1.243[N2/HZ] 
Kic=91.3  [Mpa  4m] 


L=  0.0254[  m ] 
M=  54.28  [Kg] 
Sy=560[Mpa] 

th  = 0.00254  [m] 

a=0.1354E-8 

P=2.25 


Computer  Solution  Method: 

The  first  dynamic  analysis  and  former  fatigue  fractures  analysis  are  carried  out  on  the 
basis  of  Monte-Carlo  method.  Finite  element  method  is  used  as  numerical  solution,  which 
is  carried  out  by  “NISAII”  software.  At  the  beginning,  Monte-Carlo  analysis  method  is 
presented  for  random  solution. 

Monte-Carlo  Simulation  Method: 


This  method  is  the  simulation  of  an  experiment  with  a computer.  The  set  of  random 
numbers  are  generated  for  random  parameters  at  the  beginning,  then  random  numbers  are 
constituted  in  response  equations  from  which  the  set  of  random  numbers  for 
identification  of  random  behavior  is  obtained.  This  set  is  analyzed  by  statistical  methods. 
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nonetheless  it  is  possible  to  obtain  quantitative  and  qualitative  responses.  Applications  of 
this  method  are  enormous.  In  this  paper  flow  chart  of  method  is  shown  in  Fig.  3. 


Fig. 3 flow  chart  of  reliability  analysis  by  Monte-Carlo  method 


Dynamic  Analysis: 

Random  vibration  analysis  is  carried  out  on  the  basis  of  given  dynamic  data  and  white 

noise  excitation  r:  = 1 24326  [— — ] for  obtained  RMS  of  response  that  required  in 

v e ‘ HZ 

“Damage  linear  accumulation  method”  and  determination  of  excitation  frequency  band 
that  required  in  “Calculation  of  excitation  standard  deviation  for  transient  dynamic 
analysis”.  Therefore,  stress  response  Syy  is  evaluated  for  node  on  the  crack  tip.  Power 
Spectral  Density  (PSD)  of  Syy  stress  response  is  shown  in  Fig.  4 A RMS  value  for 
loading  process  (force  excitation)  and  stress  dynamic  response  is  obtained. 

(RMS)  ex.  -1.6[N]  (RMS)Rc,  = \ .5E2[Pa] 


Fig. 4 PSD  Response  of  stress  Syy  in  frequency  domain  related  to  node  on  crack  tip 
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Results  of  Crack  Growth  and  Reliability: 

Results  of  reliability  analysis  are  presented  in  Tab. I for  data  of  final  crack  length  at  time  t=l. 6 
sec,  and  calculated  reliability  value  is  written  only  for  other  times. 


Tab!  Results  of  Reliability  and  Fracture  Analysis  (ac=  5.08[mm]) 


NUMBER  OF 
SIMULATION 

1 

2 

3 

B 

5 

6 

■ 

8 

9 

10 

2-f  [mm]  at  Time®  1.6 

3.26 

3.6 

3.92 

2.76 

ES 

3.48 

ES 

3.08 

FAILURE(a=ac) 

NO 

NO 

NO 

NO 

NO 

NO 

NO  FAILURE 

10 

7 

4 

0. 

0. 

TIME(S) 

t=l  .6 

t=2.4 

t=2.8 

t=3.6 

t=4 

R(t) 

10/10 

=1 

0.7 

■ 

0 

0 

Conclusion: 

Estimated  reliability  is  shown  in  Fig.  2 by  three  methods  “Stochastic  Averaging  (SA)”, 
“Damage  Linear  Accumulation  (DLA)”  and  “Software  (NISAII)”.  Proceeding  from  Fig. 
2 one  can  see  that  the  results  of  this  method  are  agreeable  and  very  compatible  with  two 
methods  mentioned  previously.  The  new  proposed  method  enables  us  to  analyze  wider 
range  of  specimens  in  more  complicated  conditions. 
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